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Abstract 

The impact that information diffusion has on epidemic spreading has recently attracted 
much attention. As a disease begins to spread in the population, information about the 
disease is transmitted to others, which in turn has an effect on the spread of disease. In this 
paper, using empirical results of the propagation of H7N9 and information about the dis¬ 
ease, we clearly show that the spreading dynamics of the two-types of processes influence 
each other. We build a mathematical model in which both types of spreading dynamics 
are described using the SIS process in order to illustrate the influence of information dif¬ 
fusion on epidemic spreading. Both the simulation results and the pairwise analysis reveal 
that information diffusion can increase the threshold of an epidemic outbreak, decrease the 
final fraction of infected individuals and significantly decrease the rate at which the epi¬ 
demic propagates. Additionally, we find that the multi-outbreak phenomena of epidemic 
spreading, along with the impact of information diffusion, is consistent with the empirical 
results. These findings highlight the requirement to maintain social awareness of diseases 
even when the epidemics seem to be under control in order to prevent a subsequent out¬ 
break. These results may shed light on the in-depth understanding of the interplay between 
the dynamics of epidemic spreading and information diffusion. 
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1 Introduction 


Understanding how diseases spread among individuals has been a reeent area of intense inves¬ 
tigation yj]. Epidemic spreading is generally considered to be a dynamic process in which 
a disease is transmitted from one individual to another through contacts in peer-to-peer net¬ 
works. To date, a vast amount of research has focused on understanding the epidemic spread¬ 
ing phenomenon, which can be categorized into the following types: (1) raidemic spreading 
dynamics in various network structures [2], such as a scale-free network [3, a small-world 
network [5, ^ and an interdependent network |7, (2) propagation mechanisms that describe 

the epidemic ^reading process, such as the Susceptible-Infected-Recovered (SIR) model for 
influenza [ 3, 


disease 


], the Susceptible-Infected-Susceptible (SIS) model for sexually transmitted 
and the Susceptible-Exposed-Infected-Recovered (SEIR) model for rabies 


liillJ]; (3) data-driven modeling approaches that tackle the epidemic transmission 1115] by 


analyzing the available rea’ 


individual interactions [[18 


datasets, such as the sea 


1911 , and contact patterns [l20l. 12111 . 


ing laws in human mobility 




m 


The majority of the aforementioned studies focused on epidemic spreading independently, 
ignoring the fact that information diffusion of the diseases themselves also has a significant im¬ 
pact on epidemic outbreaks. Eor example, the outbreak of a contagious disease in a population 
leads to the spread of information, either through the media or friends, regarding the disease. 
This information impacts the protective measures that others may take such as staying at home, 
wearing face masks, and taking vaccinations 11221] . These preventive behavioral responses upon 
receiving information regarding the disease may reduce the size of the epidemic outbreak. This 
is supported by research that has shown that people’s behavioral responses contribute to the 
control of disease spreading, an example of which that is evident in the severe acute respiratory 
syndrome (SARS) in China in 2003 11231]. Therefore, there has been an increased interest in ex¬ 
amining the interaction between epidemic spreading and information diffusion. Seminal study 
was conducted in order to understand how the diffusion of awareness, or disease information, 
affects the spread of diseases. This was done by building a model in which the spread of both 
the disease and information in a host population was accounted for 112414271] . The results indi¬ 
cated that the interaction between these two different processes of spreading decreased the size 
of the epidemic outbreak in a well-mixed population 11241] . In some cases, enough behavioral 
changes would occur in response to the diffusion of disease information that the growing epi¬ 
demic would vanish completely, this occurred even in cases where the epidemic transmission 
rate was higher than the classical threshold 1128l. 1291]. Additionally, the interplay between in¬ 
formation diffusion and epidemic spreading was elucidated on multiplex networks, where each 
type of spreading dynamics occurs on its own layer (information diffusion on communication 
layer versus epidemic spreading on physical layer) 113014321] . The epidemic threshold as it relates 
to the physical contact layer can be increased by enhancing the diffusion rate of information 
on the communication layer. The effect that behavioral changes have on a population can be 
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explained as being affeeted by the follo wing three eharaeteristies 11261] : (i) the disease state of 
the individuals, sueh as vaeeination OSL 13411: (ii) the epidemie transmission rate and recovery 
rate BO. 1351] : (iii) the contact network structure, such as the adaptive process [l36l-l38l]. Mass 
communication scholars share similar views on the causal linkages between the spread of epi¬ 
demics and the diffusion of information regarding the epidemics. For example, the outbreak of 
severe epidemics usually attracts heavy media coverage, subsequently resulting in the follow¬ 
ing three types of responses from the public: cognitive responses, such as the attention to the 
information and an increased awareness of the problem llSOn . affective responses, such as anxi¬ 
ety, fear, or even panic 115 ill , and behavioral responses, such as the adoption of new practices in 
order to replace undesirable habits 11521] . However, these are more assumptions than empirical 
facts because it is difficult to find relevant data for such a clear-cut process and, even when the 
data are available, it is difficult to separate the unique effects of information on the control of 
epidemics from other confounding factors, such as changes in the diseases, the seasons and the 
medical treatments. 

However, due to the difficulty in separating the unique effects of information on the control 
of epidemics from other confounding factors, such as changes in the diseases, the seasons, and 
the medical treatments, and the lack of relevant data, these can be regarded as assumptions as 
opposed to empirical facts. 

Recent studies regarding the interplay between information diffusion and epidemic spread¬ 
ing have focused on the suppression of epidemic outbreaks by information diffusion. The 
occurrence of a disease prompts the sharing of information regarding the disease, which leads 
to preventive measures that inhibit epidemic spreading 112413211 . Studies have also indicated 
that when an outbreak is under control people are less vigilant in sharing relevant information 
which leads to a decrease in preventative measures and may result in a recurrence of pandemic 
epidemics. Unfortunately, evidence suggests that the second outbreak of epidemics is typically 
more deadly. For example, the spread of SARS was shown to be alleviated in early March 
2003, however, later that month showed a rapid increase, as indicated in the evolution curve of 
the probable cases of SARS (Fig. 2 in Ref. ll23ll l. This rapid increase may be directly related 
to the time in which people’s attention became more focused on the Iraq War in 2003. This is 
the central issue that is addressed in the current study. First, we demonstrate a similar outbreak 
pattern using data on the spread of avian influenza A (H7N9) in China [l39l - l4lll along with 
the diffusion of disease information. Then a mathematical model was built that describes both 
types of spreading dynamics (i.e., disease and information) as an SIS process. Results using 
the model show that information diffusion can significantly inhibit epidemic spreading. Finally, 
simulation results exhibit multi-outbreak phenomena in the epidemic spreading accompanied 
by the impact of information diffusion, which is consistent with the empirical results. Our 
findings highlight the need for the maintenance of disease awareness even during times when 
epidemics appear to be under control. 
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2 Results 


2.1 Empirical Analysis 


The information diffusion during the H7N9 outbreak in China in 2013-14 was taken to illus¬ 
trate the influence that information diffusion has on epidemic outbreaks. The epidemic spread¬ 
ing data on H7N9 infection cases was obtained from the Chinese Center for Disease Control 
and Prevention. Messages that referred to ”H7N9” or ’’avian influenza A” on Sina Weibo, 
which was the largest micro-blogging system in China {http://www.weibo.com/), were used to 
measure the existence of information about H7N9. The assumption was made that posted or re¬ 
posted a message about H7N9 indicated an awareness of the existence of the disease; otherwise 
the user was considered to be unaware of it. 

The spreading process of both the epidemic and information of H7N9 is illustrated in Fig.[TJ 
The blue circles (J) and pink diamonds (Info) represent the epidemic spreading and informa¬ 
tion diffusion. The similarity between the evolutionary trends of the two domains is obvious. 
When the epidemic broke out in Apr. 2013 and Feb. 2014 (Fig.fT]), evidence suggests that many 
people were discussing ”H7N9” or ’’avian influenza A” in weibo.com. This suggests that the 
influence of information diffusion on epidemic spreading could be quite significant. Actually, 
public responses to H7N9, including actions such as staying at home or wearing face masks, 
can affect the spread of the epidemic. Two epidemic outbreaks occurred during the period un¬ 
der consideration. Interestingly, the size of the first epidemic outbreak (Apr. 2013) was smaller 
than the second (Feb. 2014), which inversely correlated with the information domain, in which 
the number of individuals discussing the disease during the first outbreak was much greater 
than the number of individuals discussing the disease during the second outbreak. This implies 
that an increased awareness of H7N9 decreased the size of the epidemic outbreak. Research on 
news diffusion offers an alternative explanation, which suggest that as an epidemic progresses, 
the uncertainty surrounding it declines and the diffusion of the relevant information also decays 
115311 . In the current case, it is possible that there was a higher degree of uncertainty among the 
Chinese people regarding H7N9 when they first heard about it in April 2013, however, they 
had become accustomed, and thus less anxious and responsive, to H7N9 in February 2014, 
even when there were more reported cases of infected people. 


2.2 Model Description 

In order to investigate the effect that information diffusion has on an epidemic outbreak, we 
propose a network model for the interaction between epidemic spreading and information dif¬ 
fusion. In this study, we consider two states of disease spreading: susceptible (S) and infected 
(I), and two states of information diffusion: aware (-r) and unaware (-). Therefore, each indi- 


5 



vidual in the system ean be eategorized as being in one of the four states: (i) S-: the suseeptible 
individual who is unaware of the existenee of the disease; (ii) S+: the susceptible individual 
who is aware of the existence of the disease; (iii) the infected individual who is unaware 
of the existence of the disease; (iv) /+: the infected individual who is aware of the existence of 
the disease. 

Using the SIS model, the transformation among these states is illustrated in Fig.[2l The two 
types of spreading processes can be described as follows: 


• At the initial step, an individual is randomly selected and is assigned to the state /+, 
which is considered as the seed of both epidemic spreading and information diffusion. 
All other individuals are assigned to the S- state. 


• Epidemic spreading: During each time step, the infected individuals (/+ and /_) are 
capable of spreading the epidemic to their susceptible neighbors (S'+ and S-) with the 
corresponding transmission probabilities and the infected individuals (/+ and /_) could 
recover to the susceptible state with the corresponding recovery probabilities. 


/5 

ai(3 

7 

£7 


the probability that S- is infected via the /_ neighbor 
the probability that 5'+ is infected via the /_ neighbor 
the probability that is infected via the J+ neighbor 
the probability that 5'+ is infected via the J+ neighbor 
the probability that J_ recover to S- 
the probability that recover to 5*+ 


{S.I. /_/_) 

(^+/_ ^ /+/_) 
(^_/+ ^ /_/+) 
^ J+/+) 


• Information diffusion: During each time step, the individuals with an awareness of the 
existence of the disease (/+ and 5'+) have the capability to transmit the information to 
their unaware neighbors (/„ and S'.) with the probability a. Additionally, the /+ and S'+ 
individuals may become unaware of the existence of the disease with the probabilities of 
A and 5\ respectively. 


a information transmission rate 
A information fading rate {S+ —)■ S'.) 

5\ information fading rate (/+ ^ /) 


The model assumes that when a susceptible individual is made aware of the existence of the 
disease (S+), they will take protective measures to avoid becoming infected and as is denoted 
as the reduction in the probability of infection {as < 1). Individuals in the J+ state will reduce 
contact with their susceptible neighbors in order to prevent the epidemic from spreading further, 
resulting in cj/ < 1. With the assumption of the independent effect of the infected probability, 
then asi = cr^cr/ when the /+ state individuals infect the S+ state individuals. When an J+ 
individual is aware of the presence of a disease, the person will increase her/his recovery rate by 
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taking medicine or other positive measures, which is represented by factor e > 1. Additionally, 
the /+ state individuals, which could be assumed to have a better understanding of the disease, 
would be less likely to forget relevant disease information, leading to (5 < 1. In the current 
study, since the spreading processes of information and disease are primarily determined by 
the corresponding transmission probabilities, we focused on the effect of the two parameters a 
and f3, while fixing the other parameters. The other parameters are set as as = 0.3, cr/ = 0.6, 
6 = 0.8, e = 1.5, A = 0.15 and 7 = 0.1, unless otherwise noted in the following analysis. 


2.3 Model Analysis 

The proposed model is performed on a random network with a total population N = 10000 and 
an average degree (k) = 15. We denote the infected level (/) and the informed level (Info) 
as the fraction of infected individuals (both J+ and J_) and the fraction of individuals who are 
informed of the existence of the disease (both 5+ and /+). The results of the simulation of 
the epidemic spreading process under the influence of the information diffusion when /3 = 0.3 
is shown in Fig. [3l In this model, the parameter a can be considered as the informed level 
in the system, where a large a indicates that information spreads much easier resulting in an 
increase in the number of informed individuals (as in the inset of Fig. [3]). Fig. [3] shows that 
the increase in a leads to the decreased epidemic spreading rate and the overall diminished epi¬ 
demic outbreak size. Therefore, increasing the number of informed individuals and improving 
self-protection measures may be an effective strategy to inhibit the spread of epidemics, which 
is consistent with the empirical analysis shown in Fig.[TJ 

The model described in Fig. [2] indicates that there is a mutual feedback between infor¬ 
mation diffusion and epidemic spreading. A higher prevalence of infected individuals results 
in the maintenance of more informed individuals for a smaller information fading probability 
(S < 1). This leads to the high informed level in the system, which in turn inhibits the epidemic 
spreading (a^j^s.si} < !)• This feedback effect can be clearly illustrated in the full set of dif¬ 
ferential equations based on the classical mean-field analysis (See Eq.© in Sec. Method 
and Materials). Additionally, a full set of differential equations based on the pairwise analysis 
was obtained ll42l - l44l] that describes the interaction between the two spreading processes (See 
Eq.© in Sec. Method and Materials). Actually, the pairwise analysis is also based on the 
mean-field assumption for calculating the change of each type of nodes as well as the node¬ 
pairs. The simulation results of the infected dynamics, numerical results of classical mean-field 
analysis and pairwise analysis when a = 0.6 and (3 = 0.3 is shown in Fig.lH The results of the 
pairwise model are clearly a better fit to the simulation than the results of classical mean-field 
analysis. 

In order to investigate the effect of the mutual interaction between a and 13 on the spreading 
process, we explored the full phase diagram showing the fraction of infected individuals in the 
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whole population as a function of parameters a and {3 in Fig. [5l The results of the pairwise 
analysis and the simulation, which are highly consistent with each other, are shown in Fig.5 (a) 
and Fig.5 (b). The dashed curve in each plot displays the critical point in the epidemic spreading 
process, (/Sc, etc), at which above this point an epidemic outbreak will occur in the population. 
The results clearly show that the spread of disease is much more rapid under conditions with a 
larger (3 and smaller a, which also indicates that the disease information diffusion can impede 
the the spread of disease. It should be noted that the process degenerates into a standard SIS 
model when a = 0, the condition at which there is no information diffusion in the system, 

and the outbreak threshold value of the epidemic spreading is (3c = = 0.0067 iQ]. This 

{k) 

is also consistent with the results of the pairwise analysis and the simulation shown in Fig. [51 
A detailed view of the pairwise analysis (a, (3 e [0, 0.05] in Fig.5 (c)) was plotted in order to 
illustrate the threshold changes. The threshold value of (3 is about 0.0067 when a < 0.01, 
as the disease information does not spread out in this case according to the inset of Fig. [3l 
When a > 0.01, the epidemic threshold can be significantly increased for the outbreak of the 
information diffusion. In other words, the information diffusion can effectively increase the 
epidemic threshold. 

The informed level is only slightly raised with the increase in a when a is large enough 

(e.g., a > 0.3), as shown in Fig.[3l which leads to an increase in the epidemic threshold as the 

increase in the a value is not so obvious. This results also indicates that information spreading 

is not always an effective auxiliary measure in the control of the spread of disease, for example, 

in the case in which there exists a disease with a strong infectiveness (e.g., large epidemic 

transmission probability (3 in the red range in Fig. [S]), enhancing the awareness of the disease 

alone is not enough to control the spread of disease. In order to gain a better understanding 

of the critical phenomenon, we investigate the evolution of the infection densities for various 

values of /3 when a = 0.6, as shown in Fig.[^ From the differential equations of the standard 
dl 

SIS model, — = —I + (3{k)I(l — I) (where / = /_ + /+), we obtain I cc at the critical 
point with integration, which indicates that the infection density has a power-law decay along 
the time evolution at the critical point. The inset of Fig. shows a power-law decay of the 
infection density when (3 ~ 0.0444, whereas the infected density (/) tends to be a steady-state 
value that leads to an endemic state when (3 = 0.05, and rapidly decays to zero leading to 
a healthy state when (3 = 0.04. Therefore, it can be inferred that the critical value of (3 is 
approximately 0.0444 in this case, which is consistent with the results shown in Fig. [51 where 
the phase diagram indicates that the critical value of (3 is around 0.045 when a = 0.6. 


Interestingly, the empirica 


analysis illustrates that the dynamics of many diseases exhibit 


a multi-outbreak phenomena [[l4 [45l-l47[l. in which there are several waves in the process of 


epidemic spreading, similar to the dynamics illustrated in Fig. [H Generally, there are many 
complex factors that contribute to multi-outbreaks, including seasonal influence, climatic vari¬ 
ation, and incubation period. The periodic outbreaks in the SIS model can be interpreted by the 
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influence of information diffusion. As previously mentioned, there is a mutual feedback be¬ 
tween information diffusion and epidemic spreading in the proposed model. On the one hand, 
a larger proportion of infected individuals should result in an increase in preventive behavioral 
responses [l 28 n due to the increased awareness of the disease, thus leading to a steady decrease 
in infected cases over time. On the other hand, when the spread of the disease appears to be un¬ 
der control (e.g., the size of infected population decreases), people become less vigilant, which 
leads to a decrease in the dissemination of information and increases the chances of a second 
outbreak. Notably, the size of the second outbreak is often larger than the previous one, as in 
the case of SARS in 2003 (Fig.2 in Ref. 1^ 1. Dengue Fever in Taiwan in 2001-2002 (Fig.3A 
in Ref. kzD)- 

In order to illustrate the multi-outbreak phenomena of epidemic spreading with the influ¬ 
ence of information diffusion, two critical infected levels were set in this model for simplicity 
(Ihigh and Iiow)- When the fraction of infected individuals in the population is larger than the 
critical value of the high infected level (> high), the information spreads more quickly and 
there is an increase in preventative behavioral responses. Based on this assumption, we obtain 
simulation results of the dynamic of infected individuals, as shown in Fig. |71 In this case, we 
set (3 = 0.18, which is much larger than the epidemic threshold without considering informa¬ 
tion diffusion. The disease spreads very quickly at first while there are very few people who 
are aware of the occurrence of the disease. Typically, the diffusion of information regarding a 
disease occurs at a faster rate than epidemic spreading. Therefore, as the information regarding 
the disease quickly spreads, the high informed level has a significant effect on inhibiting the 
spread of the epidemic (the decay period of the epidemic). It should be noted that the epidemic 
will be completely suppressed if the high informed level can be maintained. However, when 
the disease is controlled from the first outbreak (i.e., the infected density is smaller than Iiow), 
members of the population are likely to no longer consider the disease as a threat, thus ignoring 
the disease propagation and no longer actively engaging in protective measures, which will in 
turn lead to another outbreak of the epidemic, with a high probability of the second being larger 
than the first. In this model, such changes in behavioral responses can be illustrated by using 
different parameter settings, such as small a and large a/ and as- This allows for the disease 
to spread again as long as social awareness becomes low enough, as is shown in the second 
epidemic outbreak illustrated in Fig. |71 Similar to the empirical analysis shown in Fig. [H the 
size of the first epidemic outbreak is smaller than that of the second one, whereas the informed 
level in the first epidemic outbreak is higher than that in the second one. It should be noted 
that, due to the difficulty in precisely quantifying the informed level in the empirical analysis, 
the number of tweets that discuss the related disease is used as a proxy measure in Fig. [TJ In 
contrast to the trend shown in Fig. [H the high informed level must be maintained during the 
period when the infected proportion decreases as shown by the pink curve in Fig. |71 Based on 
the model analysis, it could be concluded that it is important to strengthen public awareness 
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of disease occurrenee, espeeially during times in whieh the spread of the epidemic is under 
control, otherwise, there is a high probability of a second outbreak. 


3 Conclusions & Discussion 


In this study, we have studied the interaction between epidemic spreading and relevant infor¬ 
mation diffusion. The empirical analysis shows that information diffusion can significantly 
inhibit epidemic spreading, in which the size of the epidemic outbreak is influenced by the 
informed level. In line with previous works on the closed feed-back loop ’’epidemic spreading 
—)■ behavior change —)■ epidemic spreading” [l24l. 13411. we build a model in which the two types 
of spreading dynamics are described as an SIS process. Both the simulation results and the 
pairwise analysis reveal that information diffusion can increase the epidemic outbreak thresh¬ 
old, diminish the final fractions of infected individuals and significantly slow down the rate of 
propagation. More importantly, we address the issue of multi-outbreak phenomena of epidemic 
spreading with information diffusion as the governing. The results of the simulation suggest 
that a higher epidemic prevalence impels people to increase in disease information sharing, 
leading to a high level of informed individuals, which in turn results in the steady decrease in 
the number of infected cases. During the periods in which the disease appears to be controlled, 
less attention is given to the disease leading to a decrease in the transmission of information 
and an increase in the chance of another outbreak. Additionally, the simulation results of the 
multi-outbreak phenomena were consistent with the empirical analysis. 

The findings from this work support the idea that preventive behavioral measures brought 
about by disease information can significantly inhibit the spread of an epidemic, and the in¬ 
crease in information diffusion can be utilized as an auxiliary measure to efficiently control 
epidemics. The government should make an effort to maintain social awareness of the disease, 
even during times in which the epidemic seems to be under control, in order to prevent another 
outbreak. In this study, we focus on the inhibition of information diffusion, and preventive 
behavioral responses are illustrated with some parameters generally. However, the dynamics 
of an epidemic may be very different depending on the behavioral responses of people, such as 
adaptive process Obll . migration [l48n . vaccination ll34ll . and immunity I149n . This work provides 
a basic understanding of the interplay between the two spreading processes. Future research 
should focus on the in-depth study of preventive behavioral responses induced by diffusion of 
disease information. 
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4 Methods and Materials 


In this study, [*] represents the number of state variables (*) in the system at time step t. 
[5”+], [S-], [/+] and [/_] represent the number of aware suseeptible, unaware suseeptible, aware 
infeeted and unaware infected. In the pairwise analysis, [*] also denotes the number of the 
corresponding state variable of the edges, for example, [*5+1+] represents the number of edges 
between two individuals at states and /+. 

Mean-field Analysis: According to Fig. |2l we adopt mean-field analysis for the spread of 
epidemic and information in a homogeneous network as follows: 


^=- - {kWiPiU]- 

d[S+] 


N 


- {fcM[5+] + + A[5+] + 7[/-] 


dt 

d[I-] 

dt 

d[I+] 


= - {k)asl3[I-]^-^ - [/+]%! + (fc)a([,5+] + [/+])%! - A[^+] + e^{I+] 


N 


=(fc)/3[J_ 


[ 5 ] 




N 


(fc)cr//?[/+]-^ - (fc)a([5+] 


= {k)(Tsld[I- 


[^+] 


{k)asail3[I+] 




N 


.<5A[/+]-7[/-] 


( 1 ) 


+ (fc)a([5+] + [/+])1^ - S\[I+] - ej[U] 


dt ' ' J TV ' ' ' TV 

Where N is the number of individuals in the system, (k) is the average degree of the network 
and the other parameters are illustrated in Nomenclature. 

Pairwise Analysis: Pairwise models have recently been widely used to illustrate the dynamic 
process of ^idemics on networks, as those models take into account the edges of the net¬ 
works H-Q. In this study, we consider a set of evolution equations which are comprised 
of four types of individuals and 10 types of edges. Using the well-known closure, expressed 


as [ABC] = 


[ABJBC] 

\B] 


(assuming the neighbors of each individual obey Poisson distribu- 
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tion) Mm, 

' rf[g-] 

dt 

d[S+] 

dt 

d[I-] 

dt 

d[I+] 

dt 

d[S-I-] 

dt 


d[S-I+] 

dt 


d[S+I-] 

dt 


d[S+I+] 

dt 


d[I-I-] 

dt 

d[I-I+] 

dt 


d[I+I+] 

dt 


d[S-S-] 

dt 

d[S-S+] 

dt 


d[S+S+] 

dt 


we can get a set of differential equations as follows : 

- /3[^_/_] - ai/3[S.I+] - a([^-5+] + [5_J+]) + A[5+] + 7[/-] 

- as(3[S+I.] - asai/3[S+I+] + a([5_^+] + [S.I+]) - A[^+] + £ 7 !/+] 


:/3[5_/_] + aiP[S.I+] - ai[S+I.] + [/_/+]) + SX[I+] - j[I.] 


■■asl3[S+I.] + asaiP[S+I+] + a([5+/_] + [/_/+]) - SX[I+] - e^[I+] 






[ 5 _] 






[/-] 


= -aiP[S-I+] + A[5+/+] - a[S-I+] - (5A[5_/+] - ^ 


+ asaiP 


[^_] [5+] [^_] 

[5+/+][5_5+] , ^ [5_/_]([/_/+] + [5+/_]) ^ [^_/+]([5_/+] + [S.S+]) 


[5+] 
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-£7[^-/+]+7[f-f+] 


= -asP[S+I.] + SX[S+I+] - X[S+I.] - a[S+I.] - + /3 


+ cr//3 
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[^_/+][5_5+] , ^ [5_/_]([^_5+] + [S.I+]) ^ [5+/_]([/_/+] + [S+I.]) 
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[I-] 
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[S-] 
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[5+] 
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[I-] 


= 2/3[5'_/_] + 2SX[I-I+]] + 2/3%^ + 2aj;3 - 2a - 27 [/_/_] 
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[ 5 _] 


[I-] 
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Figure 1: (Color online) Empirieal analysis of the epidemie spreading (blue eireles) and the 
information diffusion (pink diamonds) of H7N9. 
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Figure 2: (Color online) Diagram illustration of epidemie spreading as well as information 
diffusion. 



Figure 3: (Color online) Dynamies of the epidemie spreading proeess with different a. The 
inset shows the informed level as a funetion of a. 
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Figure 4: (Color online) The mathematieal analysis of the epidemic spreading process: simula¬ 
tion (pink circle), pairwise analysis (green solid curve) and classical mean-field analysis (blue 
dashed curve). 
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Figure 5: (Color online) Comparison between pairwise analysis and simulation for the fraction 
of infected individuals in the stationary state (colors represent the fraction of infected individ¬ 
uals). (a) pairwise analysis; (b) simulation; (c) detailed view of pairwise analysis. 
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Figure 6: (Color online) Infeetion density as a function of (3 with the pairwise analysis. The in¬ 
set is the infection density as a function of time with the pairwise analysis around the threshold 
(Different colors correspond to different (3). 



Figure 7: (Color online) Multi-outbreak phenomena of the epidemic spreading with the influ¬ 
ence of the information diffusion. I^igh and Iio^ are set as 0.05 and 0.0003 respectively. 
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